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ABSTRACT 

The  Undique  imaging  Monte  Carlo  simulator  (Un- 
dique  hereafter)  was  developed  to  reproduce  the 
behavior  of  3D  imaging  devices.  This  paper  de¬ 
scribes  its  high  fidelity  imaging  algorithm. 

1.  INTRODUCTION 

The  field  of  3D  imaging  is  evolving  rapidly.  Ad¬ 
vance  in  Flash  ladar  [1-2],  simplified  the  task  of 
laser  3D  imaging.  These  technologies  can  be  used 
for  observation  of  targets  through  foliage  [1],  aer¬ 
osol  [3],  or  netting  [2]. 

We  developed  Undique  in  an  effort  to  under¬ 
stand  the  influence  of  these  phenomena  on  Flash 
lidar  image  quality.  The  simulator  reproduces  not 
only  the  scattering  of  the  photon  in  the  propaga¬ 
tion  range  but  also  3D  solid  targets  and  basic  re¬ 
ceiver  behavior. 

In  section  2  we  give  a  brief  description  of  Un¬ 
dique  capabilities.  Section  3  describes  the  high  fi¬ 
delity  imaging  algorithm.  Finally,  section  4 
presents  some  examples  of  simulated  results. 

2.  UNDIQUE 

Undique  ( un-di-kwe )  is  a  latin  word  meaning 
“from  all  parts”,  a  name  quit  adequate  when 
propagating  photons  in  aerosols.  It  is  a  versatile 
tool  developed  to  reproduce  the  behaviour  of  dif¬ 
ferent  types  of  active  imaging  system  such  as  clas¬ 
sical  lidar  system,  passive  2D  imager,  active  im¬ 
ager,  or  range  gated  3D  imager.  To  simulate  these 
systems,  we  need  a  2D  gated  detection  system,  a 
source,  and  a  propagation  range.  In  the  case  of  3D 
range  gated  imager,  a  solid  3D  target  is  also  simu¬ 
lated.  Fig.  1  shows  a  schematic  view  of  the  mul¬ 
tiple  components  emulated  by  the  system. 

2.1  Simulated  sub  systems  characteristics 

Undique  generates  individual  photons.  The  source 
is  monochromatic  but  its  wavelength  can  be  ad¬ 
justed  before  simulations.  The  photons  polariza¬ 


tion  state  is  defined  using  the  Stokes  parameters. 
To  simulate  a  gated  camera,  we  keep  a  record  of 
the  total  distance  travelled  by  the  photons. 

The  source  can  be  located  anywhere  in  the  sim¬ 
ulated  environment.  Multiple  sources  or  pre¬ 
defined  source  patterns  such  as  bar  patterns  of  ra¬ 
dial  patterns  can  be  used. 


Fig.  1.  The  Undique  Monte  Carlo  simulation  en¬ 
vironment.  We  show  the  various  components  of 
the  system. 

The  gated  camera  is  made  of  a  perfect  thin  lens 
and  a  square  detector.  The  aperture  diameter,  the 
focal  length  and  the  focus  of  the  camera  can  be 
changed.  To  provide  polarization  analysis  capabil¬ 
ity,  a  polarizer  and  a  retarder  are  located  in  front 
of  the  detector. 

The  propagation  range  is  the  distance  between 
the  camera  and  the  solid  target.  It  is  divided  in  a 
series  of  propagation  planes.  Each  plane  generates 
a  specific  behavior.  The  simulation  basic  planes 
are  the  camera’s  aperture  plane,  the  solid  target 
upper  limit,  and  the  target  background  plane.  We 
also  have  planes  to  define  the  obscurant  limits  and 
for  turbulence  effects  (still  under  development). 
Photons  propagate  from  one  plane  to  the  next. 
When  it  arrives  at  a  plane,  actions  are  taken  and 
propagation  is  modified  accordingly. 
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The  solid  target  is  at  the  end  of  the  propagation 
range.  It  is  a  3D  structure  that  rests  on  an  infinite 
solid  background.  We  can  have  up  to  20  individu¬ 
al  targets.  These  targets  can  be  cylindrical,  rectan¬ 
gular,  or  have  specific  regular  patterns.  Albedo 
and  BRDF  of  solid  targets  can  be  changed. 

Undique  relies  on  phase  function  information  to 
scatter  the  photons.  The  first  aerosols  included  in 
the  simulator  were  spherical  water  droplets  de¬ 
scribed  by  the  Mie  theory.  The  aerosols  can  be 
monodispersed  or  have  a  normal,  log  normal,  bi¬ 
log  normal,  or  gamma  distribution. 

The  simulator  can  also  use  Rayleigh  scattering, 
Henyey  Greenstein  phase  function  and  user 
defined  phase  function  for  atmospheric  propaga¬ 
tion.  In  water,  Fournier  Forand  and  pure  water 
phase  functions  are  available. 

Aerosols  can  be  located  anywhere  in  the 
propagation  range.  It  can  span  from  the  camera  to 
the  target  or  be  concentrated  in  a  determined  posi¬ 
tion.  We  can  simulate  up  to  three  clouds  having 
either  uniform,  linearly  changing,  or  parabolic  ex¬ 
tinction  profile. 

2.2  Direct  propagation  and  variance  reduction 

When  used  in  direct  propagation  mode,  the  simu¬ 
lator  propagates  ballistic  photons  until  they  are 
captured  by  the  camera  or  reach  a  predetermined 
delay.  In  order  to  increase  the  throughput  of  the 
simulator,  the  variance  reduction  mode  generates 
events  every  time  an  interaction  occurs  within  the 
system  Field  Of  View  (FOV).  Three  different  vari¬ 
ance  reduction  methods  are  used.  The  first  method 
is  the  well-known  semi-analytical  method  [4]  also 
known  as  the  expected  value  method  [5].  The 
second  is  an  importance  sampling  method  used  to 
increase  the  amount  of  events  occurring  in  front  of 
the  receiver.  The  third  method  is  splitting  and 
Russian  roulette  used  to  produce  a  defocusing  ef¬ 
fect. 

3.  HIGH-FIDELITY  IMAGE  GENERATION 

Conventional  atmospheric  Monte  Carlo  [6]  are 
usually  coaxial  with  the  source,  focused  at  infin¬ 
ity,  and  observe  symmetrical  system.  The  camera 
aperture  is  not  considered  in  the  imaging  process. 
The  errors  caused  by  these  approximations  are 
usually  less  important  than  the  intrinsic  variability 
of  the  simulation. 


In  recent  years,  active  systems  have  been  integ¬ 
rated  in  various  industrial  processes  [7].  The  aper¬ 
ture  diameter  can  represent  a  few  percent  of  the 
total  range.  In  these  conditions,  the  photon  distri¬ 
bution  over  the  camera’s  aperture  can  have  a  sig¬ 
nificant  impact.  When  simulating  engulfed  sys¬ 
tem,  the  aperture’s  contribution  is  paramount. 


Fig.  2.  Phase  function  projection  over  a  small 
(left)  and  a  large  (right)  solid  angle  Q. 


In  direct  propagation  mode,  ballistic  photons 
touching  the  camera  aperture  are  focused  on  the 
detector  respecting  the  optical  properties  of  thin 
lenses.  In  variance  reduction  mode,  a  sub  photon 
having  the  probability  to  reach  the  receiver  is  send 
toward  it.  This  probability  is  provided  by  the 
phase  function.  Fig.2  shows  the  two  possible 
cases.  On  the  left,  the  range  is  large  and  the  cam¬ 
era’s  aperture  sustains  a  small  solid  angle  Q.  On 
the  right  the  solid  angle  sustained  by  the  aperture 
is  large  and  the  projection  of  the  phase  function 
changes  considerably  over  it.  When  the  camera’s 
aperture  is  small  compared  to  the  distance  to 
travel,  the  probability  to  go  toward  the  small  solid 
angle  Q  is 
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where  F(6,  (p)  is  the  phase  function  as  a  function 
of  elevation  6  and  azimuth  (p,  d  is  the  range  and  r 
the  aperture  radius.  When  the  camera’s  aperture  is 
larger,  the  phase  function  projection  over  the  aper¬ 
ture  varies.  The  scattered  light  can  go  toward  a 
larger  Q  and  the  new  probability  is 
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This  calculation  can  be  perform  faster  by  put¬ 
ting  F(0,  (p)sinOdO  in  a  look-up-table  (LUT). 
We  can  either  produce  an  exact  mapping  of  the 
aperture’s  probability  distribution  to  locate  the 
photon  or  use  an  alternate,  less  computer  time 
consuming  strategy  We  chose  the  latter  and  we 
explain  how  we  proceeded  in  the  next  section. 

3.1  The  photon  pointing  into  the  receiver 

In  Fig.  3,  an  interaction  occurs  at  point  A  inside 
the  camera  FOV,  and  the  photon  travels  along  line 
AB;  point  B  is  inside  the  aperture.  The  circumfer¬ 
ence  of  the  aperture  is  divided  in  24  sections  of 
15°.  They  are  used  to  define  24  slices  having  point 
B  as  vertex.  Knowing  the  interaction  coordinate  in 
A,  the  contact  point  coordinate  in  B,  and  the  edge 
coordinate  of  a  slice  in  E,  we  can  find  the  eleva¬ 
tion  angle  0  defined  by  these  three  points.  Know¬ 
ing  the  azimuthal  angle  dtp  sustained  by  the  slice, 
we  can  find  the  probability  to  be  scattered  in  that 
slice  using  Equ.  2.  The  operation  is  repeated  for 
all  the  slices  and  we  get  a  probability  distribution. 
A  random  number  is  drawn  to  determine  the  slice 
touched  by  the  photon.  Using  the  LUT  as  a  prob¬ 
ability  distribution  along  elevation,  a  second  ran¬ 
dom  number  is  used  to  find  the  photon  location  in 
6  and  a  third  is  used  to  determine  tp. 


Fig.  3.  Schematic  view  of  the  photon  distribution 
on  the  camera's  aperture  in  variance  reduction: 
the  photon  pointing  inside  case.  The  slices  are 
shown  in  yellow  and  blue. 

3.2  The  photon  pointing  outside  the  receiver 

This  situation  is  shown  in  Fig.  4.  The  interaction 
between  the  photon  and  the  aerosol  happens  at 
point  A  but  this  time  the  projected  location  of  the 
photon  is  outside  the  aperture  (point  B).  The  aper¬ 


ture  is  divided  in  24  slices  having  the  same  angu¬ 
lar  width  dtp.  A  line  is  traced  between  B  and  E,  the 
far  edge  of  the  aperture.  The  probability  to  be 
scattered  between  point  D  and  E  is  computed  for 
every  slice  using  0i,  02,  and  dtp  in  Equ.  2.  From 
that  point,  the  distribution  of  the  photons  on  the 
aperture  is  similar  to  the  previous  case. 


Fig.  4.  Schematic  view  of  the  photon  distribution 
on  the  camera's  aperture  in  variance  reduction: 
the  photon  pointing  outside  the  aperture  case.  The 
slices  are  shown  in  yellow  and  blue. 

4.  SOME  EXAMPLES 

In  this  section  we  present  simulations  made  us¬ 
ing  Undique  to  demonstrate  its  versatility.  The 
first  example  in  Fig.  5  represents  a  solid  target 
made  of  concentric  cylinders  seen  through  aero¬ 
sol.  This  was  done  to  assess  the  impact  of  aerosol 
on  our  capacity  to  detect  3D  shapes  through  ob¬ 
scuration  [8]. 


Fig.  5.  Simulation  of  a  solid  target  seen  through 
an  aerosol  made  of  10  jum  water  droplets.  Top 


left:  ODO.  Top  right:  OD1.  Bottom  left:  OD3. 
Bottom  right:  OD5  (from  [8] ). 

In  Fig.  6  we  look  at  a  radial  source  through  an 
F/2  system.  The  source-camera  distance  is  equal 
to  100  times  the  aperture  diameter  and  the  FOV  is 
50  mrad.  On  the  left,  the  source  is  focused;  on  the 
right,  defocus  generates  the  expected  phase  inver¬ 
sion  in  the  radial  pattern.  This  is  only  possible  if 
the  photons  are  distributed  adequately  on  the  aper¬ 
ture. 

In  Fig.  7  we  show  the  case  of  a  polarised  signal 
propagating  through  a  test  case.  The  case  is 
chosen  because  it  can  be  calculated  analytically 
and  because  it  is  highly  sensitive  to  small  errors  in 
polarization  management. 


Fig.  6.  Phase  inversion  of  a  radial  pattern 
generated  by  defocus  of  the  optical  system. 


Fig.  7.  Polarized  source  reflected  by  aerosols  in  a 
test  case.  We  see  the  total  reflected  signal  (top 
left),  the  principal  polarizationftop  right),  the  45  ° 


polarization  (bottom  left)  and  the  secondary  po¬ 
larization  (bottom  right). 

5.  CONCLUSION 

The  Undique  Monte  Carlo  simulator  is  a  versatile 
tool  that  can  be  used  to  reproduce  most  passive 
and  active  imager  observing  a  scene  under 
obscuration.  Its  unique  imaging  algorithm  enables 
it  to  reproduce  the  behavior  of  actual  systems 
under  a  wide  range  of  conditions. 
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